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ABSTRACT 


A  unidirectional  fiberglass  epoxy  composite  plate  specimen  is  mod¬ 
elled  as  a  homogeneous  transversely  isotropic  continuum  plate  medium. 
Acousto-ultrasonic  non-contact  input-output  characterization  by  trac¬ 
ing  SH  waves  in  the  continuum  is  studied  theoretically  with  a  trans¬ 
mitting  and  a  receiving  transducer  located  on  the  same  face  of  the 
plate . 

The  single  reflection  problem  at  a  stress -free  plane  boundary  in  a 
semi- inf inite  transversely  isotropic  medium  whose  isotropic  plane  is 
parallel  to  the  plane  boundary  is  analyzed  first.  It  is  found  that  an 
incident  SH  wave  results  in  a  reflected  SH  wave  only;  the  amplitude 
ratio  of  the  reflected  SH  wave  to  the  incident  SH  wave  is  negative 
one;  and  the  angle  of  reflection  of  the  reflected  SH  wave  is  equal  to 
the  angle  of  incidence  of  the  incident  SH  wave.  The  balance  in  energy 
flux  normal  to  the  plane  boundary  is  also  checked. 

The  delay  time  is  calculated  as  if  the  SH  waves  were  propagating 
in  an  infinite  half  space.  It  is  found  that  the  directional  depen¬ 
dence  of  the  phase  velocity  of  the  SH  waves  travelling  in  the  trans¬ 
versely  isotropic  medium  has  a  significant  effect  on  the  delay  time, 
as  opposed  to  the  directional  independence  of  the  phase  velocity  of 
the  SH  wave  travelling  in  an  isotropic  medium. 
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The  displacement  associated  with  the  SH  wave  in  the  plate  detected 
by  the  non-contact  receiving  transducer  is  approximated  by  an  asymp¬ 
totic  solution  for  an  infinite  transversely  isotropic  medium  subjected 
to  a  harmonic  point  load.  The  polar  diagrams  for  the  directivity 
function  of  the  shear  stress  due  to  SH  waves  in  the  plate  are  shown  at 
frequencies  of  0.75,  1.50  and  2.25  MHz.  To  maximize  signal  detection, 
the  transverse  transmitting  transducer  should  be  placed  such  that  the 
applied  force  generated  by  it  is  in  the  direction  perpendicular  to  the 
line  connecting  the  transmitting  and  the  receiving  transducers.  This 
study  enhances  the  quantitative  understanding  of  acousto-ultrasonic 
nondestructive  evaluation  (NDE)  parameters  such  as  the  stress  wave 
factor  (SWF)  and  wave  propagation  in  fiber  reinforced  composites  or 
any  other  materials  which  can  be  modelled  as  transversely  isotropic 
media. 
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INTRODUCTION 


Fiber  reinforced  composite  materials  are  attractive  materials 
for  aerospace  applications  because  of  their  high  specific  mechanical 
properties.  It  has  been  shown  that  many  composites,  such  as  the  uni¬ 
directional  fiberglass  epoxy  composites  or  fiber  reinforced  ceramics, 
as  shown  in  Fig.  1,  may  be  modelled  as  a  homogeneous  transversely  iso¬ 
tropic  continuum  [1] .  In  this  work,  acousto-ultrasonic  (AU)  non- 
contact  input-output  characterization  of  a  homogeneous  transversely 
isotropic  elastic  plate  is  investigated  by  tracing  SH  waves. 

First,  single  reflection  problem  of  an  incident  SH  wave  at  a 
stress-free  plane  boundary  in  a  semi- inf inite  transversely  isotropic 
medium  whose  isotropic  plane  is  parallel  to  the  plane  boundary  is  con¬ 
sidered.  At  such  boundaries,  the  conditions  for  the  existence  of  wave 
mode  conversion,  the  angle  of  reflection  of  the  reflected  wave,  and 
the  amplitude  ratio  of  the  reflected  wave  to  the  incident  wave  are 
derived. 

Second,  the  SH  wave  input -output  relations  are  derived  when  mul¬ 
tiple  reflections  occur  at  the  top  and  bottom  faces  of  the  plate.  The 
delay  time  between  input  and  output  versus  the  distance  separating  the 
transmitting  and  receiving  transducers  is  analyzed.  The  directivity 
function  of  the  shear  stress  associated  with  the  SH  waves  is  computed. 
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And,  the  output  displacement  at  the  non-contact  receiving  transducer 
is  approximated  by  an  asymptotic  solution. 

This  investigation  should  enhance  the  quantitative  understanding 
of  AU  NDE  parameter  such  as  the  stress  wave  factor.  It  also  provides 
the  potential  for  assisting  in  the  development  of  better  NDE  schemes 


utilizing  the  SWF. 


SINGLE  REFLECTION  PROBLEM  AT  A  STRESS -FREE  PLANE  BOUNDARY  IN  SEMI- 
INFINITE  TRANSVERSELY  ISOTROPIC  MEDIUM  WHOSE  ISOTROPIC  PLANE  IS  PARAL¬ 
LEL  TO  PLANE  BOUNDARY  FOR  INCIDENT  SH  WAVE 


1 .  Reflected  SH  Wave 


For  a  homogeneous  linearly  elastic  transversely  isotropic  con¬ 
tinuum,  the  number  of  independent  elastic  constants  is  five  [1] . 

Define  a  coordinate  system  (x,  y,  z)  for  a  semi- inf inite  transversely 
isotropic  medium  whose  isotropic  plane  is  parallel  to  the  plane  bound¬ 
ary  where  the  reflection  occurs  as  follows:  the  plane  boundary  con¬ 
tains  the  X  and  y  axes ,  and  the  z  axis  is  the  zonal  axis  of  the 
medium,  which  is  in  the  direction  parallel  to  the  fiber  direction  shown 
in  Fig.  1,  as  shown  in  Fig.  2.  The  generalized  Hooke's  law  is  writ¬ 
ten,  relative  to  the  (x,  y,  z)  coordinate  system,  as  [1] : 

’’xx  ”  Ciiu,x  +  C2_2V,y  +  Ci3W,z 
’■yy  ”  +  C3]^v,y  +  C]^3W,z 

’’zz  -  Ci3U,x  +  Ci3V,y  +  C33W,z 

rxz  "  C44(u,z  +  w,x)  (1) 

^yz  ”  C44(v,z  +  w,y) 

’’xy  “  C66<u.y  +  v,x) 
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where  r^g  (r,  s  -  x,  y  and  z)  are  the  normal  (r  -  s)  and  shear  (r  ^  s) 
stresses  with  respect  to  the  coordinate  system  (x,  y,  z) ;  u,  v  and  w 
are  the  displacement  components  of  a  point  in  the  medium  along  the  x, 
y  and  z  axes,  respectively;  denotes  partial  differentiation  with 

respect  to  the  variable  which  follows;  and  C-^2>  *^13 >  *^33  >  ^44 

are  the  five  independent  elastic  constants  where  Cgg  —  1/2(0]^]^  -  C-^2)  • 


Let  a  plane  progressive  wave  be  represented  as  [2] 

(u,  V,  w)  -  A  (P^,  Py,  Pz)  expCiuCSj^x  +  Syy  +  S^z  -  t) )  (2) 

where  Sj^,  Sy  and  are  the  components  of  the  slowness  vector,  which 
is  in  the  same  direction  as  the  normal  to  the  wavefront  and  whose  mag¬ 
nitude  is  equal  to  the  reciprocal  of  the  magnitude  of  the  phase  velo¬ 
city  [1],  along  the  x,y  and  z  axes,  respectively;  Pj^,  Py  and  P^  are 
the  components  of  a  unit  vector  of  particle  displacement  along  the  x, 
y,  and  z  axes,  respectively;  A  is  the  amplitude  of  particle  displace¬ 
ment;  t  denotes  time  and  w  denotes  radian  frequency.  It  follows  from 
Eqs.  (1)  and  (2)  that  the  stresses  can  be  represented  as 
r^x  “  i‘^A[CiiSxPx  +  Ci2SyPy  +  C13S2PZ ]  exp{  iwCS^x  +  Syy  +  S^z  -  t) ) 

Tyy  “  iwA[C2^2Sx^x  ^ll^y^y  Cj^^S^Pz]  ®xp { iw(Sj^x  +  Syy  +  S^z  -  t) ) 

r^z  -  iwA[Ci3SxPx  +  Ci3SyPy  +  C33S2Pz]exp{ iwCS^x  +  Syy  +  S2Z  -  t) ) 

Tjjz  -  iwA[C44S2Px  +  C44SxPz]exp{iw(SxX  +  Syy  +  S2Z  -  t) }  (3) 

Ty2  ”  f^A[C44S2Py  +  C44SyP2 ] Gxp { iw( SjjX  +  Syy  +  S^z  -  t) ) 

xy  “  ^^^[^66^y^x  ^  CggSjjPy] exp{  iw(Sj^x  +  Syy  +  S^z  -  t) ) 


The  stress  boundary  conditions  on  the  stress -free  plane  boundary 
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require  that  [2] 


r  (I)  +  ,  (R)  .  0 

^J2.  ^  ^yz  ^ 


^  ,„(R)  -  0 


where  (r-x,  y  and  z)  represents  stresses  on  the  plane  boundary 

associated  with  the  incident  wave,  and  (r  -  x,  y  and  z) 

represents  stresses  on  the  plane  boundary  associated  with  the 
reflected  waves.  In  order  to  satisfy  Eq.  (4),  it  is  required  [2]  that 
the  frequency,  w,  of  the  reflected  waves  be  equal  to  that  of  the 
incident  wave  and  that 


S,(I)  -  s,(R) 
Sytl)  -  Sy<R) 


As  a  result  of  Eq.  (5),  the  slowness  vectors  of  the  incident  and 
reflected  waves  lie  in  a  plane  called  the  plane  of  incidence.  This 
analysis  can  be  simplified  by  assuming  that  the  slowness  vectors  of 
the  incident  and  reflected  waves  are  in  the  plane  x  -  0,  as  shown  in 
Fig.  2.  Since  the  slowness  vectors  of  the  incident  and  reflected 
waves  are  in  the  plane  x  -  0, 


Sx^^)  -  -  0. 
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It  has  been  shown  that  an  SH  wave  with  slowness  vector  in  a 
plane  containing  the  zonal  axis  of  a  transversely  isotropic  medium 
possesses  a  transverse  displacement;  that  is,  for  the  coordinates  in 
Fig.  2,  (Px,  Py.  Pz)  ”  (1-  follows  from  Eqs .  (3)  and 

(6)  that  the  stresses  associated  with  SH  waves  are 

’’xx  ’  ^yy  “  ’’zz  “  ’’yz  ”  ° 

^xz  ^  0  :  ^xy  ^  0- 

It  has  also  been  shown  [3]  that  an  SV  wave  and  a  P  wave  with  slowness 
vector  in  a  plane  containing  the  zonal  axis  of  the  transversely  iso¬ 
tropic  medium  are  quasi -transverse  and  quas i- longitudinal ,  respec¬ 
tively;  the  component  of  the  unit  vector  of  particle  displacement  for 
SV  and  P  waves  along  the  x  axis,  P^,  vanishes,  whereas  the  components 
along  the  y  and  z  axes,  Py  and  P^.  do  not.  It  follows  from  Eqs.  (3) 
and  (6)  that  stresses  associated  with  P  and  SV  waves  are 

’’xy  “  ’’xz  -  0  (3) 

’■yz  °  !  ’’xx  ''  0  :  ^yy  ^  ^  ’  ’’zz 

Assume  that  an  SH  wave  is  incident  on  the  plane  boundary ,  the 
x-y  plane  in  Fig.  2.  It  follows  from  Eqs.  (4)  and  (7)  that 
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r  (I) 

xz 


'zz<^>  -  0 


'xz»>  '  0. 


As  a  result  of  Eq.(9),  it  is  known  from  Eq.  (4)  that  and 

are  equal  to  zero.  This  means  that  no  P  or  SV  waves  will  be 
reflected  back  into  the  medium  because  a  reflected  wave  of  either  the 
P  or  SV  type  results  in  nonzero  values  of  the  stresses  Tyy  and  ’"zz- 
follows  from  Eq.  (4)  that 


r  (I)  +  r  <R)  -  0 

^xz  ^  ^xz  ^ 


Since  a  reflected  SH  wave  results  in  nonzero  values  of  shear 
stresses,  and  r^y,  it  is  therefore  concluded  that  only  an  SH  wave 

will  be  reflected  back  into  the  medium. 


2 .  Slowness  Surface  for  SH  Wave 


The  equations  of  motion  relative  to  the  coordinate  system 
(x,y,z)  are  [2] 

^xx,x  ’^xy.y  ’^xz.z  “  P^>tt 

^xy,x  ^yy,y  ^  ^yz.z  “ 


’’xz.x  ^yz,y  ^zz,z  ”  P^.tt 
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It  follows  from  Eqs .  (1),  (2)  and  (11)  that  the  following  equa¬ 


tions  of  motion  are  obtained; 

t*^ll^x^  ‘  ^*^12  ^66^^x^y^y 

+  (Ci3  +  C44)SxS2Pz  -  0 

(Ci2  +C66)SxSyPx  +  [CggSx^  +  CiiSy2  0448^2  -  p]Py  (12) 

+  (C12  +  C44)SyS2Pz  -  0 

(^13  ^44)^x^z^x  ^*^13  ^44^^y^z^y  ’*’  t^44(^x  "*■  ^y  ) 

+  €335^2  _  0 

The  condition  of  the  plane  wave  solution  is  expressed  by  setting 
the  determinant  of  the  matrix  of  the  coefficients  of  P^^,  Py  and  ?2 
Eq.  (12)  equal  to  zero  [1]: 

CCllSx^  +  C66Sy2  +  C44S,2  -p]  (C^j  ^  C66)SxSy  (C^  +  044)5^8, 

<^12  +  C46)S^Sy  ^  c^^s,2  ♦  +  C44)SyS^  =  ( 

<Ci3  +  C44)S^S,  (Ci3  +  C44)SySj  ^  jyZ)  +  0338,2  .^] 

(13) 

Expanding  Eq.  (13),  we  obtain  three  sheets  of  slowness  surface.  The 
slowness  surface  for  an  SH  wave  is  [3] 

C66(Sx2  +  Sy2)  +  0448^2  .  ^  ,  (14) 
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In  the  present  study,  the  slowness  vectors  of  the  incident  and 
reflected  SH  waves  are  confined  to  the  y-z  plane,  as  shown  in  Fig.  2. 
The  intersection  of  the  slowness  surface  with  the  plane  x— 0  is  given 
by  [3]  as 


The  numerical  values  of  the  constants  Cgg  and  p  for  the 

fiberglass  epoxy  composite  shown  in  Fig.  1  are  given  in  [1]  as 
C44  -  4.422  X  10^  N/m2,  Cgg  -  3.243  x  10^  N/m^,  and  p  -  1850  kg/m^ . 

One  quadrant  of  the  slowness  surface  for  an  SH  wave  travelling  in  the 
fiberglass  epoxy  composite  is  obtained  by  substituting  the  numerical 
values  of  the  elastic  constants  and  the  density  into  Eq.  (15),  and  is 
shown  in  Fig.  3. 


Angle  of  Reflection 


It  follows  from  Eq.  (5)  that  the  component  of  the  slowness  vec¬ 
tor  of  an  incident  SH  wave  along  the  y  axis,  Sy^^^,  is  equal  to  that 
of  the  reflected  SH  wave,  Sy(^) .  As  result  of  Eqs.  (5)  and  (15),  the 
relationship  between  the  z -components  of  the  slowness  vectors  of  an 
incident  SH  wave  and  the  reflected  SH  wave  is 


n 


(16) 


s,(i)  -  -  s,(R). 

The  minus  sign  is  due  to  the  fact  that  the  slowness  vector  of  an 
incident  SH  wave  points  out  of  the  medium,  whereas  the  slowness  vector 
of  the  reflected  SH  wave  points  into  the  medium,  as  shown  in  Fig.  2. 
Consequently,  the  value  of  the  z-component  of  the  slowness  vector  of 
an  incident  SH  wave,  S2^^\  is  negative,  whereas  that  of  the  reflected 
wave,  S^^^^  is  positive. 

The  angle  of  reflection  is  defined  as  the  angle  between  the 
slowness  vector  of  a  reflected  SH  wave  and  the  normal  to  the  plane 
boundary  where  the  reflection  occurs.  Similarly,  the  angle  of  inci¬ 
dence  is  defined  as  the  angle  between  the  slowness  vector  of  an  inci¬ 
dent  SH  wave  and  the  normal  to  the  plane  boundary,  as  shown  in 
Fig.  2.  Therefore,  the  angle  of  reflection  of  a  reflected  SH  wave  is 

9^  -  tan-l(Sy(R)/S2(^>)  (17) 

and  the  angle  of  incidence  of  an  incident  SH  wave  is 

-  tan’^(Sy^^V-Sz^^^)  .  (18) 

It  follows  from  Eqs.  (5),  (16),  (17)  and  (18)  that  the  angle  of  inci¬ 
dence  is  equal  to  the  angle  of  reflection,  as  shown  in  Fig.  2. 
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Amplitude  Ratio  of  Reflected  Wave  to  Incident  Wave 

The  amplitude  ratio  of  the  reflected  SH  wave  to  the  incident  SH 
wave  can  be  computed  from  the  boundary  condition  oh  the  stresses, 

Eq.  (4).  Assume  an  SH  wave  of  unit  amplitude  travelling  in  the  plane 
X  =  0  is  incident  on  the  plane  boundary  at  the  origin,  as  shown  in 
Fig.  2.  The  shear  stress,  associated  with  the  incident  SH 

wave  on  the  plane  boundary  can  be  obtained  from  Eqs .  (3)  and  (6)  as 

-  i‘^C44S2<^)  .  (19) 

Similarly,  the  shear  stress,  associated  with  the 

reflected  SH  wave  on  the  plane  boundary  is 

-  iwA(R)C44S2(^).  (20) 

It  follows  from  Eqs.  (10),  (16),  (19)  and  (20)  that  the  amplitude  of 

the  reflected  SH  wave  is 

a(^)  -  -1.  (21) 

Therefore,  the  amplitude  ratio  of  the  reflected  SH  wave  to  the  inci¬ 
dent  SH  wave  is  negative  one . 
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5 .  Balance  in  Energ^v  Flux  Normal  to  Plane  Boundary 


The  balance  in  energy  flux  normal  to  the  plane  boundary  must  be 
satisfied  and  is  expressed  as  [3] 

F  (I)  +  F  -  0  (22) 

z  z 

where  F_(^^  and  F_(^)  are  the  components  of  the  energy  fluxes  of  the 

Z  Z 

incident  and  the  reflected  SH  waves  along  the  z  axis,  respectively. 

The  component  of  the  energy  flux  associated  with  an  SH  wave 
along  the  z  axis  is  [4] 

F,  -  .  (23) 

It  follows  from  Eqs .  (16),  (21)  and  (23)  that  the  sum  of  the  z-compo- 
nent  of  the  energy  flux  of  the  incident  SH  wave  and  that  of  the  re¬ 
flected  SH  wave  is  equal  to  zero.  Accordingly,  the  requirement  for 
the  balance  in  energy  flux  perpendicular  to  the  plane  boundary  where  a 
reflection  occurs,  Eq.  (22),  is  satisfied. 

The  z- components  of  the  energy  flux  of  an  incident  SH  wave  of 
unit  amplitude  and  the  reflected  SH  wave,  Eq.  (23),  in  the  semi¬ 
infinite  transversely  isotropic  medium  are  shown  in  Fig.  4,  exclusive 
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of  the  radian  frequency  w,  as  functions  of  the  angle  of  incidence  with 
the  numerical  value  of  the  elastic  constant  given  in  [1]  for  the 
unidirectional  fiberglass  epoxy  composite  shown  in  Fig.l  as  = 

4.422  X  10^  N/m^.  The  balance  in  energy  flux  normal  to  the  plane 
boundary,  Eq.  (22),  is  also  shown  in  Fig.  4,  which  is  independent  of 
frequency. 
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ACOUSTO -ULTRASONIC  NON-CONTACT  INPUT- OUTPUT  CHARACTERIZATION  OF  FIBERGLASS 


EPOXY  COMPOSITE  PLATE 


It  has  been  shown  [1]  that  the  unidirectional  fiberglass  epoxy 
composite  shown  in  Fig.  1  may  be  modelled  as  a  homogeneous  trans¬ 
versely  isotropic  continuum.  For  the  axes  shown  in  Fig.  1,  the  iso¬ 
tropic  plane  of  equivalent  continuum  lies  in  the  midplane  of  the  plate 
[1].  A  cartesian  coordinate  system  (x,  y,  z)  is  chosen  so  that  the  x-y 
plane  is  the  isotropy  plane;  thus,  the  upper  and  lower  surfaces  of  the 
plate  are  at  z-  h/2  and  z-  -h/2,  respectively,  where  h  is  the  plate 
thickness.  The  properties  of  the  equivalent  continuum  model  of  the 
unidirectional  fiberglass  epoxy  composite  plate  are  [1] 

h  -  0.1m 

Cii  -  10.581  X  10^  N/m^ 

C^2  -  4.679  X  10^  N/ra2 

Ci3  -  40.741  X  10^  N/m2  (24) 

C33  -  4.422  X  10^  N/m2 
-  3.243  X  10^  N/ra2 
p  -  1850  kg/m^ 

Non-contact  transmitting  and  receiving  transducers  are  located  on 
the  same  face  of  a  fiberglass  epoxy  composite  plate  specimen,  as  shown 
in  Fig.  5.  The  unidirectional  fiberglass  epoxy  composite  plate  spe- 
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cimen  is  considered  as  a  plate  of  thickness  h  and  of  infinite  planar 
(x-y)  extent.  The  input  electrical  voltage  to  the  transmitting  trans¬ 
ducer  is  V^(t)  and  the  output  electrical  voltage  from  the  receiving 
transducer  is  VQ(t)  where  t  represents  time.  The  transmitting  trans¬ 
ducer  converts  an  input  electrical  voltage  into  a  stress,  whereas  the 
receiving  transducer  converts  a  stress  associated  with  stress  waves 
travelling  in  the  plate  into  an  output  voltage.  In  the  following  ana¬ 
lysis,  only  the  SH  waves  are  traced.  The  SH  waves  which  are  generated 
by  the  transmitting  transducer  located  above  point  0  experience  mul¬ 
tiple  reflections  at  each  face  of  the  plate,  and  then  reach  the 
receiving  transducer  located  above  point  M,  as  shown  in  Fig.  6.  Since 
the  isotropic  plane  lies  in  the  midplane  and  is  parallel  to  both  the 
top  and  the  bottom  faces  where  multiple  reflections  occur,  the  angle 
of  reflection  of  the  reflected  SH  wave  is  equal  to  the  angle  of  inci¬ 
dence  of  an  incident  SH  wave  for  each  reflection  at  each  face  of  the 
plate.  Accordingly,  the  SH  wave  travelling  from  the  input  0  to  the 
output  M  may  be  considered  as  waves  propagating  in  a  semi- inf inite 
transversely  isotropic  medium  and  travelling  to  point  M'  as  if  there 
were  no  bottom  face,  as  shown  in  Fig.  6. 

1 .  Delay  Time  and  Phase  Velocity 

Let  the  input  0  and  the  output  M  lie  in  the  y-z  plane.  Assume 
the  number  of  reflections  at  the  bottom  face  experienced  by  the  SH 
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wave  in  travelling  from  the  input  0  to  the  output  M  is  n,  as  shown  in 
Fig.  6.  With  respect  to  the  z  axis,  the  angle  of  incidence  of  the  SH 
wave  at  each  face  of  the  plate  is  6 ,  and  the  total  distance  travelled 
by  the  wave  is  Rj^.  From  the  geometry  in  Fig.  6, 

e  -  tan'l(i/2nh)  (25) 

where  i  is  the  separation  distance  between  the  input  0  and  the  output 
M,  and 


R^-i/sintf.  (26) 

The  delay  time  tn  for  the  wave  to  reach  the  receiving  transducer  is 

tn  -  Rn/Ci(<^)  (27) 

where  C]^(5)  is  the  directionally  dependent  phase  velocity  of  the  SH 
wave,  and  which  for  the  unidirectional  fiberglass  epoxy  composite 
plate  is  [1] 

Ci((?)  -  [(Cggsin^tf  +  C44Cos2(?]V2  (28) 

where  C44,  Cgg  and  p  are  given  by  Eq.  (24). 
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The  phase  velocity  as  a  function  of  the  angle  of  incidence  6 
is  shown  in  Fig.  7.  The  delay  time  is  then  computed  when  the  number  of 
reflections  n  at  the  bottom  face  of  the  plate  is  equal  to  10,  100,  300 
or  500.  The  numerical  results  are  shown  in  Fig.  8  where  the  delay  time 
t^  is  plotted  as  the  ordinate,  and  the  dimensionless  quantity  Z/h.  is 
plotted  as  the  abscissa  for  values  of  zero  to  300. 

2 .  Displacement  Detected  By  the  Receiving  Transducer 

The  displacement  detected  by  the  non- contact  receiving  trans¬ 
ducer  above  point  M,  radiated  by  the  transverse  transmitting  trans¬ 
ducer,  is  assumed  to  be  equivalent  to  the  displacement  at  point  M' , 
associated  with  the  SH  wave  propagating  in  a  semi- infinite  trans¬ 
versely  isotropic  medium.  Except  for  the  reflection  coefficients  at 
each  face  (to  be  discussed  later),  the  displacement  is  computed  as  if 
there  were  no  bottom  boundary,  as  shown  in  Fig.  6.  The  shear  stress 
at  point  M'  is  approximated  by  the  far- field  asymptotic  solution  for 
large  of  an  infinite  transversely  isotropic  medium  subjected  to  a 
harmonic  point  load. 

Now  consider  an  infinite  transversely  isotropic  medium  in  which 
the  z  axis  of  a  rectangular  cartesian  system  0(x,y,z)  is  the  zonal 
axis  of  the  medium  and  the  x-y  plane  coincides  with  the  isotropic 
plane,  as  shown  in  Fig.  9.  The  equations  of  motion  including  the 
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body  force  are  [5] 


’’xx,x  ^xy,y  ■*"  ’’xz.z  +  pX  -  pu,tt  (29) 

^xy,x  ^yy,y  ^yz,z  +  pX  —  pv,tt  (30) 

^xz,x  ^yz,y  ^zz,z  “  pw.tt  (31) 

where  r^-g  (r,  s  -  x,  y,  z)  are  the  normal  (r  -  s)  and  shear  (r  s) 
stresses  with  respect  to  the  chosen  coordinate  system  0(x,  y,  z) ;  u, 

V  and  w  are  displacement  components  of  a  point  in  the  medium  along  the 
X,  y  and  z  axes,  respectively;  X,  Y  and  Z  are  the  components  of  the 
body  force  along  the  x,  y  and  z  axes  respectively;  p  is  the  density;  t 
is  time;  and  denotes  partial  differentiation  with  respect  to  the 
variable  which  follows. 

If  we  differentiate  Eq.  (30)  with  respect  to  x  and  Eq.  (29)  with 
respect  to  y,  we  find  upon  substitution  of  the  resulting  equations  and 
using  the  appropriate  stress-strain  relations,  Eq.  (1),  that 

A,tt  --^|^(A,xx  +  A,yy)  +-^-K,zz  +  Y,x  -  X,y  (32) 

where  A  -  v,x  -  u,y. 
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For  a  harmonic  point  load  at  the  origin,  we  take  body  forces  of 
the  form 


X  =  Xo5(x)S(y)S(z)e-i“t 
Y  -  Yo5(x)S(y)S(z)e-i‘^t 
Z  =  Zo5(x)5(y)5(z)e-i‘"t 

where  5(r)  (r  —  x,  y,  z)  are  the  Dirac  delta  functions  and 
Zq  are  the  magnitude  of  the  respective  point  body  forces, 
expressed  as  a  threefold  Fourier  integral  [5] 


(33) 


^o  *  ^o 
A  can  be 


A(x,y,z,t)  - 


ACSjj.Sy.S^.t)  exp{iu)(SxX+Syy+S2Z-t))dSxdSydS2  (34) 


where 


A  (  Sjj  ,  Sy  ,  S  J 


,t)  -  l/8»r- 


A(x,y,z,t)  expliwC-S^x-Syy-S^z+t) }dxdydz. 


Similarly,  Y,x  and  X,y  can  be  expressed  as  threefold  Fourier  integrals 


Y,x  -  iwSjj 


Y  expdwCSj^x+Syy+S^z-t)  IdS^dSydS^ 


(35) 


X,y  -  iwSy 


X  exp{  iwCSjjX+Syy+S^z-t)  jdSj^^dSydS^ 
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where 


00 


X  -  l/8n- 


Xo5(x)5(y)5(z)e-i‘^t 


exp { iu>(Sj^x+Syy+S2Z "  t)  jdS^dSyCiS^ 


Xq/Stt- 


and,  similarly, 


Y  =  Y^/Sn^;  Z  -  Zq/Sit- 


It  follows  from  Eqs .  (32),  (33),  (34)  and  (35)  that 
_  ^  i(S,.Y,  -  SyX^) 

8x^uG(Sjj ,  Sy ,  S2) 

where  G(Sjj,  Sy,  S^)  “  Cg6/p(Sjj^  +  Sy^)  +  -1. 

Therefore,  G(Sx,  Sy,  S^)  -  0  represents  the  slowness  surface  for 
an  SH  wave  by  comparison  with  Eq.  (14).  As  a  result  of  Eqs.  (34)  and 
(36) ,  A  becomes 


A(x,y,z,t) 


^^x^o'^^y^o 
8jr^«G(S^,  Sy ,  S^) 


exp{ iw(Sj^x+Syy+S2Z- 1)  jdS^^dSydS^  (37) 
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The  asjrmptotic  solution  at  large  distance  from  the  point  load  is 
obtained  by  applying  the  theory  of  residues,  the  method  of  stationary 
phase,  and  the  radiation  condition  [5]  as 


A(x,y,z,t)- 


ia;(s/Y^  -  S/X^) 
2fl-R|VG|  (  |Kq|  )  V2 


exp{  i«(Sjj*x+Sy*y+S2*z- 1)  ) 


(38) 


where  Sy*,  S^*)  are  the  points  on  the  slowness  surface  for  an  SH 

wave,  G  (Sj^,  Sy,  S^)  ”0,  where  the  normal  is  parallel  to  the  OP 
direction;  R  is  the  distance  from  the  origin  0  where  the  point  load  is 
applied  to  the  point  of  interest  P  in  the  medium;  VG  is  the  gradient 
of  G  (Sjj,  Sy,  S2)  and  |VG|  is  the  magnitude  of  the  gradient  of  G(Sjj, 
Sy,  S2)  and  is  expressed  as 


|vg|  -  [G,Sx  +  G,Sy  +  G,S2]V2 


(39) 


and  is  evaluated  at  points  (Sj^*,  Sy*,  S2*)  of  the  slowness  surface  for 
an  SH  wave,  G(Sjj,  Sy,  S2)  -  0,  where  the  normal  is  parallel  to  the  OP 


direction;  Kq  is  the  gaussian  curvature  of  the  slowness  surface  for  an 
SH  wave  and  is  expressed  as 


1 


(G?Sx  +  G^S  +  G^S^) 


2  2 
Z{G  ,  (G  ,  SjjSjj  G,  SySy-G,  S^Sy) 


+  2G,Sx  G,Sy(G,SjjS2  G,  SyS^-G ,  SjjSy  G.S^S^)) 


where  S  is  the  sum  with  respect  to  cyclic  permutation  of  S^,  Sy  and 
S2,  and  (1Kq|)^/^  is  the  square  root  of  the  magnitude  of  the  gaussian 
curvature  of  the  slowness  surface,  G(Sjj,  Sy,  S^)  “  0,  and  is  evaluated 
at  points  (S^^*,  Sy*,  S^*)  of  the  slowness  surface  where  the  normal  is 
parallel  to  the  OP  direction,  as  shown  in  Fig.  9. 


The  displacement  component  along  the  x  and  y  axes,  u  and  v,  due 
to  an  SH  wave  can  be  obtained  by  direct  integration  of  the  definition 
of  A  in  Eq.  (32),  in  conjunction  with  pages  14  and  15  of  [6]  and  is 
asymptotically  expressed  as 


u(x,y,z,t)  - 


Q  *2y  c  *c  *y 

Oy  Aq  Oy  1q 


27rR(s/2  +  S  *2)|vg1(1Kc|)^/2 


exp{iu)(SjjX  +  Syy  +  S^z 


v(x,y,z,t)  - 


c  '^q  *Y  4_  q  ^2y 
■^x  ^y  ■''^o  ^x  ^o 


25rR(Sx*2  +  S„*2)|vg|(|Kg1)V2 


*  *  * 

exp{iw(Sj^x  +  Syy  +  S^z  -  t) 
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It  is  known  from  the  definition  of  A  in  Eq.  (32)  and  Eq.  (37) 
that  SH  waves  have  no  contribution  to  the  displacement  component  along 
the  z  axis,  w. 


However,  it  should  be  noted  that  there  is  only  one  point  on  the 
slowness  surface,  G(Sy^,  Sy,  S^.)  =■  0,  where  the  normal  is  parallel  to  a 
given  direction  OP.  This  is  because  the  slowness  surface  for  an  SH 
wave  is  a  spheroid,  Eq.  (14).  Consequently,  if  (x,  y,  z)  are  the 
coordinates  of  a  given  point  in  the  medium,  there  will  be  only  one  SH 
wavefront  passing  through  it  corresponding  to  a  point  (Sj^  ,  Sy" ,  S^") 
on  the  slowness  surface  where  the  normal  is  parallel  to  the  given 
direction.  Note  that  in  this  analysis,  it  has  not  been  necessary  to 
use  Eq .  (31) . 


3 .  Directivity  Function 


The  shear  stress,  associated  with  SH  waves  reaching  the 

point  M'  in  Fig.  6  is  given  from  Eqs.  (1)  and  (41)  as 

iC44«s/(Sy*2Xo  -  s/Sy*Yo) 

r„_(x,y ,  z ,  t) - ^ - 

2;rR(s/2  +  Sy*2)  |vg|  ( 1x^1 )  V2 

exp{ia>(Sjj*x  +  Sy*y  +  S^*z  -  t) )  (42) 
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The  directivity  function  Dgjj  of  the  shear  stress, 
associated  with  the  SH  wave  is  determined  from  the  amplitude  of  the 
shear  stress  in  Eq.  (42)  by  setting  Xq  -  Yq  -  R  -  1; 


C44Sz*|Sy*2  -  S/Sy*|u, 

25r(Sx*2  +  Sy*2)|vG|(|KG|)V2 


(43) 


Since  the  point  0  and  the  point  M'  in  Fig.  6  are  located  in  y-z 
plane,  the  evaluation  of  the  directivity  function  Dgjj  is  simplified  by 
setting  Sjj  -  0  in  Eq.  (43). 

Due  to  axial  symmetry  with  respect  to  the  zonal  axis,  the  z 
axis,  of  the  transversely  isotropic  medium,  the  value  of  the  directi¬ 
vity  function  thus  obtained  holds  for  all  values  of  Sj^. 

Therefore,  it  follows  from  Eq.  (39),  (40)  and  (43)  that 

C44<^S2* 

4^C66/p(C44/p)V2 

The  polar  diagrams  for  the  directivity  function  of  the  shear 
stress,  Tjjz.  associated  with  the  SH  waves  travelling  in  the  unidirec¬ 
tional  fiberglass  epoxy  composite  are  obtained  from  Eqs .  (24),  (25), 
and  (44)  at  frequencies  of  0.75,  1.50  and  2.25  MHz,  and  are  shown  in. 
Fig.  10,  where  the  angle  of  incidence  9,  given  by  Eq.  (25),  is  used  to 
determine  the  direction  in  which  the  value  of  a  point  (0,  Sy  ,  S^  )  on 


^66  *  ,  ^44 

( -  Sy*)2  +  ( - 

I  P  P 


Sz  ) 


*n2 


(44) 
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the  slowness  surface  where  the  normal  is  parallel  to  the  given  direc¬ 
tion  is  thus  obtained. 

4 .  Assumptions  on  the  Transducers 

The  non- contact  transmitting  transducer  in  Fig.  5  is  assumed  to 
transform  an  electrical  voltage  into  a  uniform  shear  stress;  however,  the 
non-contact  receiving  transducer  in  Fig.  5  transforms  a  displacement  into 
an  electrical  voltage.  The  approach  which  follows  is  similar  to  that 
given  in  [ 7 ] . 

If  an  input  voltage  of  amplitude  V  and  frequency  w  is  applied 
according  to 

Vl(t)  -  (45) 

see  Fig.  5,  the  stress  a  that  is  introduced  into  the  specimen  by  the 
non-contact  transmitting  transducer  is 

a(t)  -  Fi(w)Ve-i("t  +  "^l)  (46) 

where  F]^(w)  is  the  transduction  ratio  for  the  transmitting  transducer 
in  transforming  a  voltage  to  a  stress  and  ((>1  is  a  phase  angle.  In 
Eqs.  (45)  and  (46),  the  harmonic  character  of  the  signals  is  expressed 
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in  the  complex  notation  where  i—  7-1  and  only  the  real  parts  of  these 
and  subsequent  equations  should  be  considered.  Thus,  the  amplitude  T 
of  the  applied  force  is  defined  as 

T  =  Fi(w)V.  (^7) 

Similarly,  if  a  stress  wave  producing  a  displacement  component  U 
and  frequency  w  that  is  detected  by  the  non- contact  receiving  trans¬ 
ducer  is  defined  as 

u(t)  -  Ue^'^^  (^8) 

the  output  voltage  from  the  non-contact  receiving  transducer,  see 
Fig.  5,  is 

Vo(t)  -  F2(w)U' *^2)  (^9) 

where  F2(«)  is  the  transduction  ratio  for  the  non- contact  receiving 
transducer  in  transforming  a  displaceraerft  to  a  voltage,  and  <f>2  is  a 
phase  angle.  Thus,  the  amplitude  V'  of  the  output  electrical  voltage 
is 

V'  -  F2(w)U.  (50) 
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The  characteristics  of  F]^(w)  and  F2(w)  are  unknown  except  that  the 
dimension  of  the  product  F]^(w)F2(w)  is  [kg/ra^  •  sec^]  . 

5 .  Steady-State  Output  Voltage  Amplitude  due  to  Multiple  Wave  Reflec¬ 
tions  in  Plate 


Since  only  the  SH  waves  travelling  in  the  y-z  plane  are  consid¬ 
ered  in  Fig.  6,  it  follows  from  Eqs.  (26),  (39),  (40)  and  (41)  that 
the  amplitude  of  the  hypothetical  displacement  at  point  M'  is  Uj.^'  ,  and 
is  defined  as 


(C44P)^/2  s/  t 

Um- -  (51) 

27r  (C]_]_-C3^2) 

where  T  represents  the  magnitude  of  the  shear  stress  generated  by  the 
non- contact  transmitting  transducer  along  the  direction  perpendicular 
to  the  line  connecting  the  point  0  and  the  point  M  in  Fig.  6  and  is 
equivalent  to  Xq  in  Eq.  (41). 

However,  with  the  bottom  boundary  present,  the  wave  is  reflected 
a  total  of  (2n-l)  times,  as  shown  in  Fig.  6.  Thus,  the  amplitude  of 
the  displacement  at  point  M  is  Uj|,  and  is  obtained  by  modifying 
Eq.  (51)  as 


2n-l 

%  ”  Qshsh  %' 


(52) 
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where  QshsH  reflection  coefficient  of  SH  waves  to  SH  waves,  and 

its  value  is  -1,  as  given  by  Eq.  (21). 

The  amplitude  of  the  output  voltage  from  the  non-contact  receiv¬ 
ing  transducer  is  V'  and  can  be  obtained  by  combining  Eqs .  (50), 
(51)and  (52)  as 


F2(w)  QshsH  (C44P)^/^  S^*  T 

V' - 

hr  Rn  (Cii-Ci^2) 

Substitution  of  Eq.  (47)  into  Eq.  (53)  gives 

2n-l 

Fi(a,)  F2(a))QsHSH  (C44P)^/^  s/  V 

V' - . 

hr  (C11-C12) 


(53) 


(54) 


Introducing  the  SH  wave  attenuation  constant  a  of  the  fiberglass 
epoxy  composite  and  a  possible  electrical  signal  amplification  factor 
K,  Eq.  (54)  can  be  written  as 

2n-l 

K  Fi(«)  F2(«)  Qshsh  (C44/’)^/^  Sj,*  V 

V' - (55) 

2ir  Rn  (C]^]^-C2^2) 

Eq.  (55)  gives  the  output  voltage  amplitude  from  the  non-contact 
receiving  transducer  due  to  an  input  voltage  amplitude  V  at  the  non- 
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contact  transmitting  transducer  when  the  SH  wave  path  has  included  n 
reflections  from  the  bottom  face  of  the  fiberglass  epoxy  composite 
plate  specimen,  as  shown  in  Fig.  6. 
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DISCUSSION  AND  CONCLUSION 


In  the  acousto -ultrasonic  input- output  characterization  of  the 
unidirectional  fiberglass  epoxy  composite  plate,  the  angle  of  reflec¬ 
tion  of  the  reflected  SH  wave  is  equal  to  the  angle  of  incidence  of 
the  incident  SH  wave  for  each  reflection  at  either  the  top  or  the  bot¬ 
tom  face  of  the  plate.  This  is  due  to  the  fact  that  the  isotropic 
plane  is  parallel  to  either  face  of  the  plate.  However,  if  the  paral¬ 
lelism  between  the  isotropic  plane  and  the  plane  boundaries  where 
reflection  occurs  does  not  exist,  the  angle  of  reflection  is  not  equal 
to  the  angle  of  incidence.  In  such  a  case,  the  use  of  a  semi- inf inite 
transversely  isotropic  medium,  neglecting  the  existence  of  the  bottom 
face  of  the  plate  (except  for  the  cumulative  reflection  coefficients) , 
to  compute  the  delay  time  t^  and  the  directivity  function  Ds^  becomes 
inappropriate . 

It  is  observed  from  Fig.  8  that  an  increase  in  the  number  of 
reflections  n  from  the  bottom  face  of  the  plate  results  in  a  minor 
increase  in  the  delay  time  t^  at  each  value  of  dimensionless  separa¬ 
tion  i/h.  For  a  given  plate  thickness  h  and  a  given  separation  dis¬ 
tance  i  between  the  transmitting  transducer  and  the  receiving  trans¬ 
ducer,  an  increase  in  the  number  of  reflections  n  results  in  a 
decrease  in  the  angle  of  incidence  8,  as  given  by  Eq.  (25),  thereby 
increasing  the  travelling  distance  R^,  from  Eq.  (26).  The  phase 
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velocity  Ci  of  an  SH  wave  in  the  unidirectional  fiberglass  epoxy 
composite  increases  with  decreasing  angle  of  incidence  6  as  shown  in 
Fig.  7.  The  increase  in  the  travelling  distance  is  partly  offset 
by  the  increase  in  the  phase  velocity  C]_  with  an  increase  in  the 
number  of  reflections  n.  This  accounts  for  the  minor  increase  in  the 
delay  time  given  by  Eq.  (27).  On  the  other  hand,  in  an  isotropic 
medium,  the  phase  velocity  of  a  stress  wave  is  directionally  indepen¬ 
dent.  Thus,  the  increase  in  the  delay  time  tj^  caused  by  an  increase 
in  the  number  of  reflections  n  is  solely  attributed  to  the  increase  in 
the  travelling  distance  in  an  isotropic  plate.  Therefore,  it  is 
concluded  that  the  directional  dependence  of  the  phase  velocity  of  a 
stress  wave  travelling  in  a  transversely  isotropic  medium  has  a  sig¬ 
nificant  effect  on  the  delay  time  when  conducting  acousto-ultrasonic 
testing. 

Now,  consider  the  case  of  a  given  number  of  reflections  n.  An 
increase  in  the  separation  distance  i  results  in  an  increase  in  the 
travelling  distance  R^,  as  shown  in  Fig.  6,  and  an  increase  in  the 
angle  of  incidence  9  given  by  Eq.  (25).  The  increasing  angle  of 
incidence  9  results  in  the  decreasing  phase  velocity  as  shown  in 
Fig.  7.  Consequently,  the  delay  time  t^j  increases  sharply  as  a  result 
of  the  increasing  separation  distance  i  for  a  given  number  of  reflec¬ 
tions  n,  as  shown  in  Fig.  8.  However,  as  the  separation  distance  2 
approaches  infinity  ,  the  phase  velocity  approaches  a  limit.  Thus, 
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the  increase  in  the  delay  time  t^^  is  due  solely  to  the  increase  in  the 
travelling  distance  R^,  which  is  similar  to  the  isotropic  medium  case. 

An  examination  of  Eqs .  (41)  shows  that  the  y-component  displace¬ 
ment  V  associated  with  an  SH  wave  travelling  in  the  y-z  plane,  ■=  0, 
is  zero  and  the  x- component  displacement  associated  with  an  SH  wave 
travelling  in  the  y-z  plane,  —  0,  is  attributed  solely  to  the 
x-component  of  the  applied  point  load.  In  other  words,  the  non- 
contact  receiving  transducer  is  unable  to  detect  the  displacement  com¬ 
ponent  along  the  line  connecting  the  transmitting  and  the  receiving 
transducers  by  tracing  SH  waves  in  the  plate.  Further,  the  only 
detectable  displacement  component  which  is  perpendicular  to  the  normal 
to  both  faces  of  the  plate  and  to  the  line  connecting  the  transmitting 
and  the  receiving  transducers  (that  is,  the  u  displacement)  becomes 
non-detectable  if  the  shear  stress  introduced  by  the  non-contact 
transmitting  transducer  is  in  the  direction  parallel  to  the  line  con¬ 
necting  the  transmitting  and  the  receiving  transducers,  regardless  of 
how  large  the  stress  is.  Accordingly,  to  maximize  the  displacement 
detected  by  the  non-contact  receiving  transducer,  the  transmitting 
transducer  should  be  placed  such  that  the  shear  stress  generated  by  it 
is  in  the  direction  perpendicular  to  the  line  connecting  the  transmit¬ 
ting  and  the  receiving  transducers  as  shown  in  Fig.  11. 
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This  theoretical  investigation  provides  a  step  forward  in  the 
quantitative  understanding  of  acousto-ultrasonic  non- destructive  eval¬ 
uation  (NDE)  parameters  such  as  the  stress  wave  factor  (SWF)  in  trans¬ 
versely  isotropic  media.  It  also  provides  the  potential  for  assisting 
in  the  development  of  more  efficient  and  more  revealing  NDE  schemes 
utilizing  wave  propagation. 
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Slowness  vector 
of  incident  SH 
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Coordinate  system  (x,y,z)  in  analysis  of  single  reflection 
problem  at  stress-free  plane  boundary  of  semi-infinite 
transversely  isotropic  medium;  z  =  0  is  plane  boundary 
where  single  reflections  occurs,  z <  0  is  free  space,  and 
X =  0  is  plane  of  incidence. 
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Component  of  slowness  vector  in  y  direction,  Sy(s/km) 
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Schematic  of  non-contact  acousto-ultrasonlc  test  configuration. 
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Pig,  7  PhasG  vGlocity  versus  angle  of  incidence  for  SH  wave  travelling  in  unidirectional 
fiberglass  epoxy  composite  plate. 
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Fig.  8  Delay  time  versus  dimensionless  separation  £/h  with  number  of  reflections  n 
from  bottom  face  of  unidirectional  fiberglass  epoxy  composite  plate  specimen 
as  parameter. 
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Fig.  10  Polar  diagram  of  directivity  function  of  shear  stress  associated 

with  SH  waves  at  frequencies  of  0.75,  1.50  and  2.25  MHz  for  positive  y-z  quadrant. 


Fig.  11  Schematic  illustrating  maximization  of  displacement 

detected  by  non-contact  receiving  transducer  by  placing 
transmitting  transducer  such  that  shear  stress  generated 
by  transmitting  transducer  is  in  direction  perpendicular 
to  line  connecting  transmitting  and  receiving  transducers 
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